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EXPERIMENTAL AND CALCULATED HISTORIES OF 
VAPORIZING FUEL DROPS 

Ey R. J. Friem, G. L. Borman, M. M. El Wakil, 
0. A. Uyehara, and P. S. Myers 


SUMMARY 


The present report contains the results of an experimental and 
theoretical investigation of the vaporization of fuel droplets in heated 
air under atmospheric pressure. 

First, an experimental investigation of the temperature and mass 
histories of single droplets was made, with emphasis on small drops 
down to 500 microns. 

Second, a comparison of experimental histories with calculated 
temperature and mass histories was made. In connection with the cal- 
culated histories, the equations of change and the associated boundary 
conditions are given in reduced form along with the various dimension- 
less parameters that would appear in a solution of these equations. 

Third, the time taken by a droplet vaporizing in high- temperature 
air to form a mixture of fuel vapor and air of combustible strength in 
the film of the droplet and at the self -ignition temperature is cal- 
culated. The concept of physical ignition delay of single droplets is 
thus presented. In addition, the vaporization of fuel droplets in a 
spray under conditions where there exists an influence of one droplet 
on another is theoretically investigated. 

The experimental and calculated results can be summarized as follows: 

(1) An appreciable portion of the drop may be vaporized dining the 
"heating-up" period; at high air temperatures as much as 50 percent or 
more of the mass may be vaporized during this period. 

(2) Calculated and experimental times to reach the same percent 
npss vaporized differed by 20 percent on an average with the computed 
values usually smaller. 

(3) For 500-micron drops, computations including the heating-up 
period gave better agreement with experimental results than did computa- 
tions neglecting the heating-up period or assuming "peeling," that is. 
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zero thermal conductivity of the liquid. This is especially true for 
low-volatility fuels. 

(4) The computed physical ignition delay of single droplets increases 
rapidly with decreased fuel volatility and decreased air teirperature and 
is relatively insensitive to variations in total pressure. 

(5) If there is Interaction between the vaporizing drops in a spray, 
cooling of the air results, approaching, as a limit, the adiabatic 
equilibrium temperature . The higher the concentration and volatility 

of the fuel, the lower the adiabatic equilibrium temperature. 


INTRODUCTION 


Fuel which is injected into a combustion chamber by a nozzle leaves 
the nozzle orifice as sheets or ligaments which eventually break down 
into drops of varying sizes. As soon as these drops are formed they 
start to vaporize because of the increase of surface -to -volume ratio. 

As the drop vaporizes it simultaneously heats up until its temperature 
approaches asymptotically a constant temperature which is determined by 
the environmental conditions. Previous investigations of fuel ignition 
are discussed in references 1 to l6. 

w 

A detailed theoretical study of the unsteady-state portion of the 
total vaporization time for single droplets was performed in a previous 
investigation (ref. l). From this study it was estimated that the larger 
drops emanating from the injector of a jet engine reach the combustion 
zone while still in the unsteady-state or heating-up period. After the 
importance of the unsteady state had been verified theoretically a 
combined theoretical and experimental investigation was undertaken to 
determine if the unsteady-state period was important (ref. 2). This 
combined investigation served to check the accuracy of the theoretical 
calculations . 

The combined theoretical and experimental investigation showed that 
the theory used for the calculations produced histories which gave 
reasonable agreement with the experimental histories . The range of 
conditions covered by this first experimental investigation was limited 
to drop sizes in the neighborhood of 2,000-micron initial diameter and 
to a maximum air temperature of 1,080° R. In order to bring the experi- 
mental conditions closer to actual jet conditions the investigation was 
extended to smaller drop sizes, ■ higher air temperatures, and different 
fuels. The results of this extended investigation are reported herein. “ 

In addition to the extension of the range of conditions covered t 

experimentally, comparisons were made with calculated histories to 
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determine and to illustrate the agreement between theoretical and exper- 
imental results. By plotting the comparisons against various parameters 
the trends of the theoretical accuracy are determined, thereby indicating 
the error that would be obtained by extrapolating the theoretical tech- 
nique to conditions other than those covered in this investigation. 

Tn order .better to orient the above study to actual spray conditions 
some thought was given to the possible range of conditions found within 
the spray. Under spray conditions the assumption of no influence of one 
drop upon another may not be entirely valid. Two extreme cases may be 
encountered. The first of these is the case of the free drop; the other 
is the condition known as adiabatic equilibrium. In the case of adiabatic 
equilibrium the droplets cool the surrounding air, a condition that may 
be possible in the initial core of the spray. Actual spray conditions 
elsewhere probably fall somewhere between these extreme conditions. 

Besides the problem of the effect of one drop upon another, there 
is the added problem that under most actual spray conditions the air 
temperature is high enough to cause self-ignition of the drops before 
they are completely vaporized. With the assumption of no influence of 
one drop upon another, a theory was developed which predicts the effect 
of various parameters upon the time required to form a combustible 
mixture. This time is defined as the physical-ignition-delay period. 

The present investigation was conducted at the University of 
Wisconsin under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. 


EXPERIMENTAL APPARATUS AND TECHNIQUES 


The experimental temperature-time curves were obtained by hanging 
a drop of fuel on the junction of a thermocouple which was connected to 
a high-speed recorder. A stream of hot air was suddenly blown past the 
stationary drop (fig. l) . Hie velocity of the air was measured by a 
velometer. Hie radius -time curves were obtained by photographing the 
drop with a movie camera (fig. 2). The mass-time curves were then 
obtained by calculating the mass from the measured radii and the known 
densities at the measured temperature s . 

To obtain the small dropB discussed in this investigation it was 
necessary to make extremely small thermocouples. This proved to be a 
very difficult task. A thermocouple made of low-thermal-conductivity 
wire with a ratio of initial drop diameter to wire diameter of 20 to 1 
or greater was used to minimize the conductive heat transfer. The 
accuracy of the temperature measurements as determined by the variation 
in the steady-state temperature was ± 4 °. 
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The problem of placing a drop on the thermocouple limited the 
investigation to drops having an initial diameter of approximately 
500 microns or more. With smaller drops there was insufficient time 
to transfer the drop to the thermocouple because it either evaporated 
completely on the wire or was too large to hang on the thermocouple. 

The extremely short lifetimes of the small drops necessitated the 
use of a high-speed, direct-recording oscillograph. The instrument 
used has a response of 0.01 second for full-scale deflection and a 
maximum chart speed of 10 inches per second. Droplet histories with a 
total lifetime of 0.1 second could be recorded with this instrument. 

The air heater (fig. 1) consists of tubular electric heaters with 
a max imum permissible sheath temperature of 2,060° R. The heat losses 
in the calming section were reduced by placing this section inside the 
heater. This location also reduced the temperature fluctuations in the 
test-section airstream. A maximum air temperature of 1,440° R with air 
flows up to 6 pounds per hour was obtained with this design. 

Hie optical system for viewing and photographing the drops is shown 
in figure 2. The light source used in previous investigations was 
replaced by a Pointolite lamp which provided a greater intensity of 
illumination. A half-silvered mirror was used to provide two images. 

One image was photographed with the l6-millimeter movie camera and the 
other image was projected on a ground glass for visual observation. The 
magnification of the optical system was increased to 50 for studying 
small drops. The photographic system was calibrated by introducing a 
wire of known diameter at the location of the drop. The error for the 
500-micron drops was estimated to be ±3 percent for the initial radius 
measurements and 5 percent for the final diameter. 


EXPERIMENTAL RESULTS 


More than 300 experimental histories have been obtained and studied. 
These histories covered 10 paraffin hydrocarbons ranging in volatility 
from n-hexane to n-octadecane. The air temperature was varied between 
680° and 1,430° R; the air velocity, between 27 and .228 inches per 
second; and the initial drop diameter, between 490 and 2,180 microns. 

Drop histories were followed to a minimum drop size of 150 microns, that 
is, a size slightly greater than the bead size of the thermocouple. 

In order to illustrate the type of information obtained figure 3 
is presented. This figure gives liquid temperature, mass, and radius 
histories for small drops of different fuels. Figure 3(a) presents 
data taken with a low air temperature. Figure 3(b) presents data taken 
with a high air temperature. 
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The radius curves for high-molecular-weight fuels first increase 
and then decrease -while the curves for low -molecular -weight fuels decrease 
immediately because the high-molecular-weight fuels initially expand 
thermally while not losing much mass by vaporization. It will also be 
noted in figure 3 that there was a variation in initial liquid temperature 
with fuel volatility. In general , no correction was made for this differ- 
ence in initial liquid temperature. 

Because of the large number of histories obtained it was decided to 
present the effect of the parameters studied, that is, air temperature, 
drop size, and fuel volatility, by cross plots of results. The points 
chosen for cross -plotting were the times required to reach the steady 
state and to vaporize 20 , 50 , and 80 percent of the original mass of the 
drop. The time to reach the steady state (i.e., the unsteady state) Is 
ihe time required for the drop to reach its wet-bulb temperature. 

Figure 4 presents cross plots for the parameter of air temperature 
against vaporization time for the fuels n-hexane, n-decane, and 
n-hexadecane . The data show that the time to reach the steady state 
decreases with higher air temperature and fuel volatility. Ihe time 
required to reach the steady state also tends to approach a constant 
value asymptotically; this trend is more evident for the fuels with 
higher volatility. 

Figure 5 presents cross plots for the parameter of fuel volatility 
as expressed by the normal boiling point at air temperatures of 900° 
and 1,430° R. The data show that all of the parameters plotted are less 
affected by fuel volatility at higher air temperatures. The trend Is 
accentuated because the same initial liquid temperatures were not obtained 
for all fuels. 

Cross plots of the parameter of initial drop size are presented in 
figure 6 for the fuels n-hexane, n-decane, and n-hexadecane. While it 
is not clearly evident in figure 6 , other data show more clearly that, 
within the experimental accuracy, the time required for 20 -, 50 -, and 
80 -percent mass transfer and the unsteady-state time plot as straight 
lines in this figure. The slope of these lines (neglecting the unsteady- 
state time) was theoretically determined (ref. 2 ) as 2.0 for high air 
velocities and 1.5 for still air. Inasmuch as the slopes of the lines 
of 20 - , 50 -, and 80 -percent mass transfer are all less than 1.5 It would 
appear that the effect of the unsteady-state time was to decrease the 
dependency on drop size. While not shown here, data similar to those 
in figure 6 but for lower air temperatures show a slightly greater 
dependency on initial drop diameter. 

It can also be seen by the change in the location of the uns teady - 
state line with fuel volatility (see figs. 6 (a) and 6 (c)) that, with a 
heavy fuel, appreciably more mass is transferred during the unsteady state. 
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The steady-state temperature of a particular fuel vas dependent on 
only the air temperature for the range of conditions covered, as illus- 
trated in figure 7* The shaded area for each fuel represents the varia- 
tion in the experimental data. The effect of size or velocity on the 
steady-state temperature vas less than the experimental accuracy. 

The times required to reach the steady-state temperature and various 
percents of the mass are enlightening and shov the importance of the 
effect of various parameters on the vaporization period, but they do not 
indicate the relative importance of the heating-up period. Consequently, 
the percent of the mass that had been vaporized at the end of the heating- 
up period vas plotted as a function of the air temperature (fig. 8). The 
shaded area represents the extremes of the variations obtained. All the 
experimental points fell vithin this area irrespective of fuel, size, 
velocity, and initial temperature; therefore, any effect of these variables 
on the percent of the mass transferred at the end of the heating-up period 
vas less than the experimental accuracy. Figure 8 shove that at high air 
temperatures 50 percent or more of the vaporization process occurs during 
the unsteady-state period. 


COMPARISON OF EXPERIMENTAL AND CALCULATED RESULTS 


Histories obtained under the condition vhich gave the minimum 
deviation betveen experimental and calculated histories are given in 
figure 9(a), vhile histories shoving the maximum deviation are given in 
figure 9(b). The histories for the large n-hexadecane droplet vaporizing 
in an airstream of 1,250° R (fig. 9(a)) shov excellent agreement betveen 
experiment and theory. The variation betveen experiment and theory for 
maximum deviations (fig. 9(b)) is vithin the experimental accuracy except 
for the temperature histories vhere the greatest deviation vas 15 ° at the 
"knee" of the theoretical curve. 

The agreement betveen calculated and experimental steady-state 
temperatures is shovn in figure 10 for three of the fuels studied. The 
difference betveen the tvo results is less than the experimental varia- 
tion, vhich vas estimated to be ±t°. 

A conparison of the time required for 20-, 50-, and 80-percent mass 
transfer as obtained experimentally and as calculated is shovn in fig- 
ure 11. The shaded areas represent the regions in vhich all the points 
fell. The experimental accuracy vas approximately ±10 percent. Any 
variation vith air temperature or velocity vas smaller than the experi- 
mental accuracy. The data indicate that at large drop sizes the experi- 
mental vaporization time vas slightly less than that calculated, •vhile 
for small drops the experimental time vas approximately 20 percent greater 
than that calculated. 
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For comparison purposes calculations were performed considering the 
drop to he at the steady-state temperature throughout its lifetime. The 
* solid lines of figures 11(a), 11 (b), and 11(c) show the ratio of experi- 

mental to calculated times for 20-, 50-, and 80 -percent mass transfer, 
respectively. For high-volatility fuels (n-hexane), almost the same 
agreement is obtained with or without consideration of the unsteady-state. 
For low- volatility fuels (n-hexadecane), a considerable error occurs when 
the heating-up period is neglected. These figures also show that the 
error in not considering the unsteady-state period increases rapidly as 
the drop size is reduced. 

The effect of fuel volatility on the ratio of experimental to cal- 
culated time required to vaporize a fixed percentage of mass is better 
shown in figure 12, where the ratios of experimental to calculated vapor- 
ization times are plotted for an initial drop size of 500 microns. In 
figure 12(a), the three lower lines indicate the small variation obtained 
between the experiment and the calculation when the unsteady-state period 
is considered in the calculation. This figure shows a large variation 
with fuel volatility and clearly indicates the need for considering the 
unsteady-state vaporization period for low -volatility fuels. 

■* A comparison of the time required to reach the steady state as 

determined by experiment and by calculation is given in figure 13 as a 
function of the normal boiling temperature of each fuel. The solid 
4 vertical lines represent the variation in the ratio of experimental time 

to calculated time. The figures clearly Indicate that for high -volatility 
fuels the theory is in considerable error. For heavy fuels the agreement 
is very good, as was indicated in the a ample history given in figure 9(a). 

There is yet another technique for calculating vaporization times. 
This technique, commonly known as peeling, assumes that none of the heat 
arriving at the liquid surface goes to heat up the liquid drop; that is, 
there is zero liquid thermal conductivity. Therefore, the heat reaching 
the liquid surface just equals the heat needed to heat the vaporizing 
mass to the vaporizing temperature and to vaporize it. This vaporizing 
or surface temperature is lower than the steady-state temperature obtained 
experimentally . 

Figure 12(b) shows the comparison of experimental results with cal- 
culations based on the peeling technique. This calculation technique 
produced better agreement than the one considering no -unsteady state as 
given in figure 12(a). However, when using this calculation technique, 
the time for 20-percent vaporization is still very sensitive to fuel 
volatility, and for low- volatility fuels the error is quite large. 

It should be remembered that in all the previous comparisons the 
air velocity was held constant during the lifetime of the drop. Similar 
comparisons would be found, however, if variable air velocities were 
considered. 
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DROPLET VAPORIZATION CALCULATIONS UNDER 
CONDITIONS FOUND IN SPRAYS 


Vaporization conditions encountered in actual fuel sprays may be 
somewhat different from the conditions experimentally studied in this 
report. Reference 1 gave a calculation technique that may be used to 
include the variable droplet velocity as predicted by the drag coeffi- 
cients and flow theory. Other considerations may be needed, however, in 
order to understand combustion in a spray better. Some thought was 
given to two of these considerations and the resulting theory was used 
to predict the possible effect of various parameters briefly. 

The first situation considered is a drop injected into an atmosphere 
hot enough to cause self-ignition where there is no influence of one drop 
on another. The second situation considered was the cooling of the air 
by the vaporizing drops. The results of computations for drops vaporizing 
under these conditions will be presented below. The symbol list is given 
in appendix A and the theory and equations used to obtain these results 
are given in appendixes B and C. This theory was originally applied to 
Diesel engines and more extensive data and comparisons are given in 
reference 3* 

Speaking of the first situation, as a droplet vaporizes in a high- 
temperature atmosphere a layer of fuel-air mixture builds up around the 
drop. The fuel-vapor-air mixture in the film varies in composition 
along radial lines from the drop, with the richest mixture being next 
to the droplet surface. The temperature in the film increases from the 
liquid temperature to the air temperature . As the droplet heats up, its 
vapor pressure increases until a combustible mixture is first found at 
the liquid surface. As the drop is heated to yet higher temperatures, 
this mixture moves away from the drop into a region of higher temperature . 

When the steady-state liquid temperature is reached the vapor film 
will be fixed in both its partial pressure and temperature distribution. 

It can therefore be seen that If a combustible mixture (say stoichio- 
metric) Is to be formed around the drop at the self-ignition tenperature 
it must occur before or at the steady state with the beginning of the 
steady state as the limiting time. 

The time necessary to form a chemically correct mixture at the self- 
ignition temperature has been defined as the physical-ignition-delay 
period. It is supposed that any further ignition delay is due to chem- 
ical factors; that is, it will be chemical ignition delay. 

The effect of air temperature upon the computed physical ignition 
delay of a drop of n-decane is shown in figure 14. The ignition delay 
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shows a sharp decrease with increasing air temperatures with a leveling 
off at higher air temperatures. The physical ignition delay goes to 
« Infinity at low temperatures, indicating that the combustible mixture, 

as previously defined, cannot form at its self -ignition temperature. 

The results of computations (not presented here) when the total 
pressure was increased from l/2 to 5 atmospheres gave only a small 
increase in physical ignition delay, indicating that the physical igni- 
tion delay is comparatively insensitive to variations in total pressure. 

The higher volatility fuels form a stoichiometric mixture more 
rapidly (as shown in fig. 15) despite their lower molecular weight and 
the resulting necessity of a higher vapor pressure to form a stoichio- 
metric mixture. It should be noted that heavier fuels reach a higher 
liquid temperature before a combustible mixture is formed and consequently 
they may experience cracking in the liquid phase before the combustible 
mixture is formed. This fact may be important in smoke and deposit 
problems . 

The second situation is concerned with adiabatic equilibrium within 
a spray. The work previously reported in this paper assumed that a large 
w excess of air always surrounded the vaporizing drop and that at no time 

did the vaporization process affect the temperature or composition of the 
air surrounding the drop except in the film. This assumption seems 
t plausible in sprays where there is considerable distance between fuel 

droplets and where the fuel-air ratio is low. If, however, the core of 
the spray Is considered it seems conceivable that the vaporization process 
may cool the surrounding air appreciably and that the fuel vapor in the 
atmosphere may reduce the driving force for mass transfer. The equations 
applicable under this assumption are given in appendix D. Briefly, It is 
assumed that there is no net heat or mass transfer to or from, the section 
under consideration and that the air and fuel vapor external to the drop 
and Its Immediate film are uniform at any time. 

Figures 1 6 and 17 show the results of calculations made using the 
theory of appendix D. Figure 16 is a plot of air -fuel-mixture tempera- 
ture versus overall fuel-air ratio. This is a thermodynamic calculation 
and thus is not concerned with the amount of time necessary to reach 
this end condition. Figure 1 6 shows that a considerable drop in tenpera- 
ture may occur even for low fuel-air ratios. Fuel-air ratios as high or 
higher than stoichiometric (approximately O. 067 ) may occur in the primary 
zone and thus it Beemed desirable to perform calculations which would 
give some indication of the time required to approach adiabatic equilibrium. 

Figure 17(a) shows histories for drops of n-heptane vaporizing with 
fuel-air ratios of zero, O.O 67 , and 1.0. The temperature of the air and 
vapor drops rapidly while the mass and liquid tenperature histories are 
L affected only slightly as compared with those of pure air when a fuel-air 
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ratio of O .067 is considered. The percent mass vaporized at any time 
decreases rapidly as the fuel-air ratio is increased. Figure 17(13) 
shows similar histories for n-hexadecane with similar conclusions. 

The above conclusions and results were based solely on the theory 
presented in appendixes C and D. Wo direct experimental evidence was 
available for comparison with the theoretical predictions. 


DISCUSSION OF RESULTS 


The results of this investigation of the experimental and calculated 
temperature and mass histories of vaporizing fuel drops have shown that 
the unsteady- state vaporization period is Important in the vaporization 
of low -volatility fuels of all size drops. For high-volatility fuels, 
such as n -hexane, where only the time for complete vaporization is of 
interest, there is less error in neglecting the heating -up period. 

The calculation techniques described in reference 1 gave reasonably 
accurate histories of the percent mass transferred from the drop. The 
discrepancy between calculated and experimental histories for large drops 
and high-molecular -weight fuels is within the experimental accuracy, but 
the discrepancy increases with decreasing drop size and fuels with lower 
molecular weights to a maximum of 20 percent for n-hexane drops of 
500-mlcron Initial diameter. This indicates that the use of the calcu- 
lation technique for very small drops might introduce a considerable 
error. The inclusion of the heating-up period in the calculations gives 
a significant Improvement over other techniques for small drops. 

The increase in error with decreasing drop size is believed to be 
due to the omission of the P ft *(t) term derived in appendix B. As the 

drop size decreases, the vaporization time also decreases, thereby pro- 
ducing a larger rate of change of pressure In P^*(t). The P^*(t) term 

would therefore increase with decreasing size and likewise the error 
introduced by neglecting it would increase. 

Although the mass histories give reasonably good agreement the tem- 
perature histories obtained with the calculation technique described in 
reference 1 do not fully coincide with the experimental histories for 
both small drops and low-molecular-weight fuels as shown by figure 9» 

Tie discrepancy with drop size is believed to be due to the omission of 
the AT-j^t) term (appendix B), which would increase with decreasing 

drop size. The discrepancy with fuel volatility could result from the 
omission of the dimensionless term where M a is the mass of the 

air and Mp is the mass of the fuel vapor. 
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If, because of their complexity, it is not feasible to perform the 
calculations including the unsteady state, it is suggested that calcula- 

* tions based on peeling be employed rather than calculations which consider 
the drop to be at the steady-state temperature throughout its lifetime. 

For example, the error in the vaporization time for 50- and 80-percent 
vaporization of a small drop of n-octadecane when using the steady-state 
calculation is l 60 and 100 percent, respectively (fig. 12(a)), while if 
the peeling calculation is used the error is reduced to 50 and 10 per- 
cent, respectively (fig. 12(b)). The error in the 20-percent-vaporization 
time is not significantly decreased by use of the peeling calculation. 

The theoretical work done on physical ignition delay indicates that 
this phenomenon may be of importance in cases where the fuel droplets are 
not completely vaporized before reaching a high-temperature zone. The 
theory again points out the importance of the unsteady- state portion of 
the histories. Drops injected in the reaction zone itself may exhibit 
behavior different from that indicated by the theory. Clearly a better 
understanding of chemical and physical delay is necessary to predict 
total ignition delay accurately. 

Adiabatic equilibrium conditions may exist in the very core of a 
\ spray and a short distance from the nozzle, but the calculations performed 

here show the effects of these conditions to be negligible unless the 
fuel-air ratio is above stoichiometric . This analysis gives another indi- 

* cation that single-droplet calculations may be applied to the major por- 
tion of the spray without great error. 


The University of Wisconsin, 

Madison, Wise., July 6, 1955* 
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APPENDIX A 
SYMBOLS 


A 0 surface area of liquid drop, sq in. 

Bo thickness of air-vapor film surrounding drop, in. 

C 1 = VfL^uA’T 

C 2 = Qv - wc pf T L 

c 3 = 

Cp specific heat at constant pressure, Btu/(lb)(°F) 

D diffusion coefficient of air-vapor system, sq in. /sec 

cL^ molar rate of diffusion of fuel vapor at liquid surface, 

f mole/(sq in.) (sec) 

d y _ diffusion velocity of component i with respect to a 

i plane moving at mass average velocity, in. /sec 

f molar mass flux vector with respect to liquid surface, 

mole/ (sq in.) (sec) 

g acceleration due to gravity, in. /sec 2 

Bjju Eusselt number for heat transfer, dimensionless 

J mechanical equivalent of heat, (in. ) (lb)/Btu 

Ji mass flux vector of component i with respect to plane 

moving at mass-average velocity, lb/(sq in.) (sec) 

K thermal conductivity, Btu/(in. ) (sec) (°F) 

Kj^ average thermal conductivity in air-vapor mixture, 

Btu/ (in. ) (sec) (°F) 

k coefficient of mass transfer, l/sec 
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k Nu 

Lq 

M 

v% 

mi 

n 

n i 

P 

P T 

Pr 

P A *(t) 

W 0 > 

Q 

%, 

Or 

Qs 

Qv O 

Qa 

R 

Re 

r 

r d 


Nusselt number for mass transfer, dimensionless 

scale factor for length, in. 
mass, lb 

mass of drop at beginning and end of increment of time, 
respectively, lb 

molecular mass of component i, lb/mole 
total number density, mole/eu in. 
number density of component i, mole/cu in. 
partial pressure, lb/sq. in. 
total pressure, lb/sq. in. 

Prandtl number, dimensionless 

change in fuel vapor pressure as a function of time 

vapor pressure gradient at liquid surface at zero time 
as function of 0, lb/cu in. 

total heat transfer from air to drop, Btu/sec 

beat received at drop surface, Btu/sec 

sensible heat received by drop, Btu/sec 

heat carried back with diffusing vapor in form of super- 
heat, Btu/sec 

heat transfer rate at zero time, Btu/sec 

heat of vaporization, Btu/sec 

universal gas constant, (in. ) (lb)/(mole) (°F) 

Reynolds number, dimensionless 
radius at any point in film, in. 
radius of drop, in. 
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Sc 

Schmidt number, dimensionless 

A 

T 

temperature in film at radius r, °R 

Sf 

t b 

air temperature at film boundary, °R 


*1 

intermediate temperature (see appendix D), °R 


T L 

temperature of liquid drop, °R 


^m 

mean tengperature in film, °R 


AT L (t) 

change in liquid temperature as a function of time, °R 

- 

V 9 > 

gradient of film temperature at liquid surface at zero 
time as a function of 0, °R/in. 


t 

time, sec 

- 

■tn 

any time increment 


V 

velocity, in. /sec 

* 

V 

reduced radial velocity, dimensionless 


v 0 

angular velocity, radian/sec 


w 

molecular weight, lb/mole 



weight of air and fuel, respectively, in atmosphere 


w 

mass vaporization rate, Ib/sec 


w 0 

mass transfer rate at zero time, lb/sec 


X 

body forces 


y 

fuel-air ratio by weight 


a 

correction factor for mass transfer, dimensionless 


a' = [Prp j (PfL 

- P fB )J log e j(P T - Pf B ) j (P T - Pfx] 

_ 

P 

dimensionless number, M a /M^ 

** 

w 



MCA TN 3988 


15 


r 

dimens ionleBs number, QvqLq^KAqTq 


e 

dimensionless number, P Q ^ p fL “ p fB 

) 


dimensionless number, ^oj (j-B ~ ^l) 


n 

P V 2 

pressure number, P — 

gP 0 


e 

angle 


~k 

heat of vaporization, Bfcu/lb 


p 

viscosity, lb/ ( in. ) (sec) 


p 

density, lb/cu in. 


Pm 

mass density of mixture, lb/cu in. 


T 

reduced surface area, dimensionless 


$ 

dissipation function, l/sec 2 


♦ 

dimensionless number, I^ w o c Po/ K o A -o 


X 

dimensionless number, P 2 /P 0 


CD 

total mass flow from drop 


V 

del operator, i-i 7 -^— + ip +■ lx 

OX^ 0 X 2 

d 

d*/dr 

reduced radial gradient, dimensionless 

5*/S0 

reduced angular gradient, dimensionless 

0 * 

reduced variable, dimensionless 


Subscripts : 



A 

component A 


a 

air 



B 


at film boundary 



fuel vapor 

liquid or at liquid surface 

scale factor 

component i 

component 1 

component 2 

total 

rate or gradient 

value with respect to space-fixed axis 


at time zero 
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APPENDIX B 


THEORY FOR DROP VAPORIZING IN HOT AIRSTREAM 


The problem of a drop vaporizing in a hot airstream covers the 
problem of simultaneous transfer of heat and mass in a flow system. 

Such systems can be accurately described by the "equations of change" 
of fluid mechanics. These equations of change are the equation of con- 
tinuity for each of the chemical species, the equation of motion, and 
the equation of energy balance. The various dimensional groups that 
should be correlated to form an experimental solution can be predicted 
from, these equations when they are written in dimensionless form. 


Mass Transfer 

The transport of mass as described in reference it- is dependent on 
(l) a gradient in the chemical potential, (2) a gradient in the total 
pressure, (3) a gradient in the temperature of the gas, and (4) trans- 
v port due to external forces. The gradients in the total pressure and 

temperature are considered to be negligibly small in regard to both 
position and time in this presentation. The transport of mass per unit 
i area because of a gradient in the chemical potential is determined from 

the following flux equation: 


Ji = n i m i d v ± 


(Bl) 


For a two-component system, assuming an ideal gas and no pressure gradient 
in the vapor tangential to the liquid surface, this equation becomes for 
a spherical drop 


f 


1 


= Ze D fd p l\ 

P 2 RT\dr / r =r 0 


(B2) 


If the total pressure, partial pressure of component 2, D, R, and T 
are all constant over the entire surface, equation (B2) can be integrated 
over the entire surface to obtain the total mass flow from the drop ctr 




r=r D 


dT 


(B3) 
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Two arbitrarily chosen scale factors (indicated by a subscript o), the 
length L 0 and pressure P D , are now introduced. It is convenient to 
introduce the following transformed (reduced) variables (indicated by 
a superscript *) in terms of these scale factors: 


p i* - p i/ p o 
r* = r/L 0 


Introducing the difference between the partial pressure of the fuel 
vapor at the liquid surface and that at infinity (or the film boundary). 


to = A, 


P T M f D / P fL “ P fB 


P2 RT 


i) L k 


'S3 

’fL “ PfR/\8r* 


dT 


(B4) 


r=r^ 


The pressure gradient 8P^/8r*, the diffusion coefficient, and the tem- 
perature are evaluated at the liquid surface. The integral in equa- 
tion (b 4) is similar to the integral which Schlichting obtained for heat 
transfer (ref. 5) • The integral in equation (B4) will be defined as the 
Nusselt number for mass transfer; that is. 





P - P 
fL fB 


/d p l*^ 
\8r* , 


dT 


r ~r r 


(B5) 


The/ functional form of this 


equation will be obtained later. 


Heat Transfer 

The transport of energy is dependent on_ (l) the temperature gradient, 

(2) the transport of thermal energy by the flux of the molecules, an d 

(3) the reciprocal process to thermal diffusion known as the "Dufour 
effect." By neglecting the Bufour effect the heat transfer to a 
spherically symmetrical liquid drop was shown in reference 2 to be 
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When the scale factor T 0 is used, T* = T/T 0 , and, introducing the 
difference in temperature at infinity (or the film boundary) and at the 
liquid surface, 



It is understood that the temperature gradient 3T*/dr* and the thermal 
conductivity K are evaluated at the liquid surface. The integral in 
equation (B8) is similar to the integral which Schlichting (ref. 5) 
called the Husselt number for heat transfer or 



(B8) 


The functional form of this equation is given in the next section. 


Air -Vapor Film 

The four equations of change needed to describe the film surrounding 
the liquid drop are: (l) Equation of continuity for component A, (2) 

equation of continuity for the mixture as a whole, (3) equation of 
motion, and (4) equation of energy balance. The equation of continuity 
for component A is given in reference 4 and the equations of continuity 
of the mixture, motion, and energy balance are found in reference 6. All 
the equations are applicable to turbulent or laminar flow. 

Equation of continuity for a mixture .- The cont inuit y equation for 
a mixture of compressible fluids is 

^ + (V • pV) = 0 (B9) 


where the operator D/Dt is the substantial derivative. 

Equation of energy balance .- The energy -balance equation for a 
perfect gas is 


■ m 


D (^ T ) . 

Dt 


1 IE + (y . KvT) + -it- <S 
J Dt gj 


(BIO) 
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Equation of motion .- For a compressible fluid, the equation of 
motion is 


P m 5X = -g7P + pj. + (V • | V(mV • ?) (Bll) 


Equation of continuity for component A .- Ihe continuity equation of 
component A for compressible fluids is 

^ = -(v • pV) - (V - jj) (B12) 


In addition to these differential equations, it is necessary to 
know the pressure and temperature dependence of the physical properties. 
These functions are denoted as follows : 


P = P (p a ,t) ' 

pp = c p( p A’ T ) 

p - p ! p a< t ) ? 
K =K(P A ,T) 

D - D(P A ,T) 


(B15) 


The following scale factors are also introduced: The length L e , 

velocity V 0 , density p Q , temperature T 0 , pressure P 0 , specific 
heat c Po , viscosity p 0 , thermal conductivity K Q , and diffusion 
coefficient D 0 . 

In the following discussion the properties Cp, p, and K are 

assumed to be constant at some appropriate mean value. With this assump 
tion it is convenient to pick the scale factors Cp Q , Po> and- Kq 

identical to these mean values so that p*, K*, and Cp* as defined 

on the following page are unity. The following reduced variables 
are thus defined: 
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Radius 

Radial velocity . . . 
Angular velocity . . . 

Density 

Time 

Pressure 

Viscosity 

Thermal conductivity . 
Specific heat .... 

Temperature 

Diffusion coefficient 

Molecular gas constant 
Radial derivative . . 


. r* = r/L 0 

. v r * = v r /v 0 

. V e * = V 0 L 0 /V Q 

p* = p n /p Q 

. t* = tV 0 /L 0 

. Pi* = Pi/P 0 

. (i* = p/p 0 B 1 

. K* = K/Kq 5 1 

°P* " = 1 

. T* = T/T 0 

D* = D/Do 

R* R Tog 

* M* “ M 72 

v o 

8* 8 

• — - ”* 

8r 8r 


The four equations of change are then rewritten in spherical 
coordinates with the further simplifying assumption that the body forces 
X s and the dissipation function 0 are negligible. 

The equation of energy balance: 


8T* 

8t* 


Ko 


yVo^cPo^oJ 


1 1 5 L 

8r V 8r*/ sin 0 86 \ 


sin 0 


8t* n 

8e 


(Blit) 


The equation of continuity: 


^*> + 


(r*) 


g[p*VKr*) 2 ] 


2 8r 


r* sin 0 80 


A [( sin e) p*V Q *j =0 (B15) 


The equations of motion: 



(Bl6) 
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St* 


M = 


gPo \ 1 SP* 


f ^0 ' 

\ 1 ( 

* 

1 

V V oVq 

/ (r*) 2 

sin 9 


- — sin 6 
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SV e *) 


S9 


S* 

Sr 


(r*) 


2 SVq* 


Se 


> + 


4- -— — ^1 — cot 9 esc 8 —(sin 9 Vq*) + 

VoW3r*[ S8 V ' 


— ^ofsin 9 V e *) + d| - ■ (r*) 2 V e *| 

sin 9 ^9 2 ' 89 Sr* e J 


(BIT) 


The continuity equation for component A: 
8*P* 1 


St* 


— |p*(r*) 2 V *1 ^-(p* Bin 8 V 0 *) + 

(r*) 2 Sr L J r* sin 9 89 v > 


/ V 

1 S* 

"(r*) 2 M*(P T *) 2 S*(Pl*/V) 

V v o L oy 

I (r*) 2 8r 

P*R*P 2 * Sr* 


r*M*(P T *) _ a*(P A */P T *) 

‘ D* 


t*r*p 2 * 


\V Q L Q Jr* sin 6 58 


89 


sin 9 


M* (Pijt*) 2 D* 8 *P a *P t * 


T*R*P 2 * Sr 


r* sin 8 


m*p t *d* S*(p a */p t *) 

Se 


t*r*p 2 * 


(bi8) 


nrhp four d im ension! ess parameters 'which appear above are related to 
well-known dimensionless numbers in the following ways: 
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( = (Jo_ LiiaJ = il 

^PoVoV Wa J ** Re 

_£s>_\ .1 

vpo^o 2 ; 11 


where 


Re 

Reynolds number 

Sc 

Schmidt number 

Rr 

Prandtl number 

n 

pressure number 


Additional dimensionless parameters might appear 'when one assumes 
special forms for the temperature and pressure dependence of the physical 
properties of equations (B13). For example, the ratio of the molecular 
weights appears if an ideal gas is assumed. 

In addition to the system of equations given above, the boundary 
conditions of the problem must be specified. The boundary conditions 
are as follows : 


For pressure: 


(a) 

p A * .0 

at 

r 

= OO 

for all times and angles 

(*>) 

Sfai - 0 

at 

r 

= 00 

for all times and angles 

(c) 

dp A* , 0 
de 

at 

r 

= r o 

for all times and angles 

(a) 

*A* ' 

at 

r 

= r O 

for an angles 



at zero time 


For temperature: 


at r = 00 for all times and angles 


(f) T* = 1 

(g) SgE-o 


at r = 00 


for all times and angles 
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(h) 

11 

0 

at 

r = r 0 

for all times and angles 

(i) 

T* = T*(t) 

at 

£ 

11 

u 

for all angles 

(j) 

T ■ ^ (e) 

at 

r = r Q 

at zero time 

For radial velocity: 




(k) 

v r * = 1 

at 

r = 00 

for all times and 0=0 

U) 

V r * = 1 

at 

¥' = 00 

for al 1 times and 0 = at 

(m) 

v r * = 0 

at 

r = r 0 

for all times and all angles 

(n) 

S*V * 

£_ = 0 

Sr 

at 

r = r 0 

for all times and all angles 

( 0 ) 

V r * = 0 

at 

all radii and times and 0 = ot/2 

For angular velocity: 




(p) 

< 

CD 

* 

II 

H 

at 

r - a> 

for all times and 0 = at/2 

(4) 

v 0 * =0 

at 

r = r 0 

for all times and all angles 

(r) 

S*v fi * 

.Sr ° 

at 

r = r Q 

for all times and all angles 

(s) 

v - 0 

at 

all radii and times and 0=0 

(t) 

£ 

* 

11 

0 

at 

all radii and times and 0 = at 


Conditions (k), (l), and (p) specify that the scale factor V G will 
he taken equal to the free airstream velocity. Condition (f) specifies 
that the scale factor T 0 will be equal to the airstream temperature. 

With these boundary conditions and the dimensionless numbers obtained 
in equations (EL4) to (Bl8), the temperature, partial pressure of com- 
ponent A, radial velocity, and angular velocity are all seen to be 
functions of: 

(1) Reynolds number 

( 2 ) Prandtl number 
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(3) Schmidt number 

(4) Pressure number 

(5) Pressure gradient at the liquid surface at zero time 

(6) Temperature gradient at the liquid surface at zero time 

( 7 ) Change in partial pressure of component A with time 

(8) Change in temperature of the liquid surface with time 

(9) Radius 

(10) Angle 

(11) Time 
That is. 


w 

3* = T(Re, Rr, Sc, n , Pj *(t), T*(t), 

PfL V (e), T Lv ( 9), r*, 9*, t*) 

(EL9) 


P* = P(Re, . . 

. t*) 

(B20) 


V r * = 7 r (Re, . 

. . t*) 

(B21) 


V 0 * = V e (Re, . 

. . t*) 

(B22) 


Since the gradients of the pressure and temperature at the liquid 
surface cannot be determined directly, it is convenient to transform 
these terms by use of equations (B4) and (Hf) : 


"U ■“•u 


P 2 RT Lq J s 


(B23) 


and 


ap a *< e > , f [ /jsgj£a\ 0 .i 

cir* lAoMfP^P o y J 



(B24) 


(B23) 


Qvq _ ~ m o 
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dT*(6) 

dr* 


f(re*) 


(B2 6) 


where r ^o^o* w Q is the mass transfer rate at zero time, and 

Qv 0 i s 'the heat transfer rate at zero time.* The dimensionless parameter 
which appears in equation (3324) can "be written 


(We\ = />WKo Y^qcjA /H yA( ma /p p \ 
Iwvj \JW WW V Vo / \ m o p t / ■ «%/ V p o/ 

RTud 

For a perfect gas = 1; therefore. 


where 


/wqR TL 0 P 2 \ 

\AoWfDP 0 p Tyl 




(B2T) 


(B28) 


In terms of these new dimensionless numbers, the functional forms 
of the temperature and partial pressure of component A become: 

T* = T(Re, Pr, Sc, Ty, P A *(t), T*(t), p, T, X, r*, 0*, t*) (B29) 

(B30) 


P* = P(Re, . . . t*) 
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Introducing the notation 


\*B - T J 



■the Nusselt numbers can be defined in the functional form as follows : 

Hjjfo -H(0, Re, Er, Sc, P A *(t), T*(t), 3, *, X, r) (531) 

k Nu = k(«. Re, Pr, Sc, n, P A *(t) , T*(t), p, X, r) (532) 


It is seen, therefore, that the Nusselt numbers are functions of 
12 parameters, 11 of which are identical in both the heat-transfer and 
mass-transfer numbers. Hie parameters \|r and X are similar to the 
parameters Z and a, given in the equations of references 1 and 2. 
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APPENDIX C 

THEORY OF SELF-IGNITION OF SINGLE DROPLETS IN HOT AIR 


According to current theories, ignition of a fuel-air mixture occurs 
in the vapor phase at a mixture ratio and temperature suitable for com- 
bustion. For a single droplet surrounded by a large amount of air this 
air-vapor mixture can be formed only in the film surrounding the liquid 
droplet. The temperature and mixture ratio at any point in the air- 
vapor film is controlled by the rates of heat and mass transfer and 
therefore by the temperature and mass histories of the droplet. 

The spherically symmetrical model given in reference 1 is used to 
calculate the physical delay time. The temperature in the air-vapor 
film increases from the liquid temperature T^ at the liquid surface 

to the air temperature T-g at some distance from the drop. The partial 
pressure of the fuel vapor decreases from P^ (the vapor pressure at 

the liquid surface) to zero at some distance from the drop. The curve 
of temperature against distance from the drop can be determined from the 
following equation for heat transfer as given in reference 1: 

Q v = iWwr 2 || - wcpf (T - Tl) 


or 


dT 

dr 


Qv + vc pf (T - T^) 

K^rtr 2 


(Cl) 


The 
from the 


curve of partial pressure against distance can be determined 
following equations and identities as given in reference 2: 
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Therefore, 



RTwP a 

DP^M^^irr^a, 


(C5) 


With the calculation technique of reference 1 it is possible to 
obtain temperature histories similar to those shown in figure 18 . Prom 
these histories and equations (Cl) and (C5 ) a relationship between the 
partial pressure of the fuel vapor and the temperature in the film is 
obtained. Illustrative partial pressure and temperature lines in the 
film, are presented in figures 18 (a) and 18 (b). 

If the partial pressure of the fuel vapor needed to form a com- 
bustible mixture is arbitrarily represented by the stoichiometric 
mixture, it is possible to determine the temperature of this mixture 
at any time. This mixture is represented by the horizontal line in 
figure 18(a). For the first curve drawn, the partial pressure of the 
liquid is below that needed to form the stoichiometric mixture. As 
the droplet heats up its vapor pressure rises and the point where the 
partial pressure in the film is such that a stoichiometric mixture is 
formed moves away from the liquid surface. 

The temperature of the stoichiometric mixture as a function of 
distance in the film can be determined from figures 18 (a) and 18 (b) . 

This temperature is plotted in figure 18(e). It should be noted that 
the temperature of the stoichiometric mixture increases with time. If 
conditions are right, the stoichiometric mixture reaches its self -Ignition 
temperature before the droplet completely vaporizes and after a chemical 
delay the drop burns. The self-ignition temperature Is represented by 
the dashed line of figure 18 (e) which intersects the curve at point 2. 

This point then determines the physical-lgnitlon-delay time as shown in 
figure 18(d). This is the earliest point in the vaporization process 
that a combustible mixture is formed at Its self-ignition temperature. 

The ignition delay 0^ is calculated by the following analytical 
technique: Combining equations (Cl) and (C5) and using from reference 2 

w - iS£ Vo p aV 


gives 

dFa _ ^m^fL^Wu^a 

dT p l(Qv " WG pf T L + wc pf T ) 


(06) 
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With the assumption that % and Cp f are constant at an average 
value throughout the film, the following substitution can be made: 


C 1 = VflX&Ui / P T 
c 2 = ^v “ VG pf T L 


C 3 = wc pf 


The values of C-^, 
film at any one instant 


C 2 , and 

of time. 


are then constants throughout the 
Equation _(C 6 ) becomes 


dP a C x dT 

pT = c 2 + c 5 t 


(CT) 


Integrating between the limits P a (r) = Pij. - Pp and Pa^ r=co ^ - p t w ^ Lere 
T(r) = T and T( r=00 ) = T B , 


T = 



Pf \C2/Cl c 2 


(C 8 ) 


This equation relates the temperature _T at some point in the air' 
vapor film to the vapor pressure Pf at the same point. This equation 

can be solved only in conjunction with the vaporization histories 
obtained by the technique described in reference 1 , since the constants 
C-l, C 2 , and are functions of the vaporization history. 

The partial pressure of any fuel corresponding to a given mixture 
strength (on a weight basis) is determined by the following equation 


Pf 




1 


(C9) 


Although it is possible to burn air-vapor mixtures over a range of 
mixture ratios, a stoichiometric mixture has been chosen to represent 
the combustible mixture. 
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A representative history of a drop of n-decane calculated "by the 
technique of reference 1 and equations (C8) and (C9) is shown in fig- 
ure 19. At point a the liquid temperature has reached the value where 
the liquid vapor pressure is just high enough to form the stoichiometric 
mixture. Fran this point the temperature of the stoichiometric mixture 
rises and approaches a constant value as the droplet approaches the 
steady-state temperature. 

Sometime along the path of the stoichiometric-mixture-temperature 
line a point is reached at which the temperature equals the self-ignition 
temperature. The self -ignition temperatures of pure hydrocarbons are 
found in the literature (among others, refs. J > 8, and 9)- The values 
obtained by different investigators vary according to the procedure, 
container, and purity of the hydrocarbon. The values given in refer- 
ence 7 'were used in this investigation although they are lower than those 
given by other investigators. The self -ignition temperature of the 
stoichiometric mixture is reached at point b in figure 19 and can be 
taken as the dividing point between the physical and chemical ignition 
delays of a single droplet. 
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APPENDIX D 

ADIABATIC SATURATION IN SPRAYS 


If fuel is injected in the form of a spray in an airstream, the 
static pressure change in the system due to vaporization cooling is very 
small. If no heat is transferred from adjacent segments in the spray 
and sufficient time is allowed, the temperature of any segment will 
ultimately reach equilibrium. Two different cases may arise depending 
upon the overall fuel-air ratio. For rich mixtures the equilibrium 
will occur in the wet region) that is, a three-phase system will result. 
In the case of leaner mixtures all of the fuel will vaporize and the 
vapor will be superheated to the equilibrium temperature. 

The wet region is defined as the condition where the equilibrium 
vapor-air ratio is less than or equal to the overall fuel-air ratio. 

For the wet region the heat balance may be written 


Q, = Qj. + Q^- 


Hie following equation can then be written for a constant-pressure 
system: 





(Dl) 


For the case where no liquid remains at the equilibrium region the 
following ' equation applies : 



CpL dT + 7 \ ± + 



(D2) 


where Tj_ is an intermediate temperature at which the liquid is assumed 
to be completely vaporized at a pressure equal to the partial pressure 
of the completely vaporized fuel in the system. Equations (Dl) and 
(D2) may be solved to obtain the final adiabatic equilibrium tempera- 
ture for any fuel-air ratio and initial conditions. 



5Z 


NACA TN 3988 


33 


While the above theory gives the conditions at equilibrium, it gives 
no information concerning what is happening in the spray as a function 
of t ime . 

In order to determine the time dependency, consider a spray of 
uniform drop size or some section of a spray where drop size is ass ume d 
uniform. Also assume that there is no net mass or heat transfer to or 
from this section. The approximation will be made that all concentra- 
tion and temperature gradients are confined to films of small but f ini te 
width surrounding each drop. Air and fuel vapor outside the film bounda- 
ries will be considered uniform throughout the section and will be 
referred to as the "atmosphere." It is further assumed that the droplets 
within the section have the same velocity relative to the atmosphere; 
this velocity may be constant or varying in time. For convenience one 
may think of each drop affecting only a volume of air determined by the 
overall air-fuel ratio (by weight). The temperature and composition of 
the atmosphere as a function of time may thus be determined, taking into 
account the effect of the heat and mass transfer on the temperature 
fuel partial pressure in the atmosphere. Extensive variables for the 
spray section may obviously be obtained by merely multiplying the single- 
drop results by the total number of drops. 

The mass-transfer equation given in reference 2 may be written as 



= S/X + _l\ 
™ \B 0 rj 



(D3) 


The pressure Ppp can then be written as 


Clearly, 


■£B " tTT 


W a^ 

Mf W D 


+ 1 



0 *) 
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Then 



As in reference 2 the mass -transfer equation may he written as 

» - Ak(Pfi ' CUT) 

where - Pf-g'j is now the reduced driving force. 

As the drop vaporizes the atmosphere is cooled and thus the average 
temperature of the atmosphere Tg must be calculated for each increment 
of time. The heat supplied by the atmosphere is the heat that enters 
the outer edge of the vapor film: 

Q = Qj, + Q s 

If the mass of air in the film is considered negligible Qg is the 
heat that goes into superheating the fuel vapor in the film. Therefore, 
for any time increment At^ 


T B(n+l) T Bn “ 


V 1 - t l) 

^a c pa + WfCpf 


n 


At n 


(D8) 


The stepwise technique described in reference 1 may then be used to solve 
the vaporization equations. The vaporization histories must of course 
predict a final temperature for the system that agrees with the final 
temperatures as calculated by the thermodynamic theory presented earlier. 


Histories using this technique are shown in figures 17(a) and 17 (b) . 
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Figure 1.- Schematic diagram of experimental apparatus 
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n-Hexadecane 6 


n-Dodecane 


n-Tetradecane 5 


n-Decane 


n-Octane 2 


n-Hexane 1 


IBjHiBW te 

IBM Ml II 


Tlme^ sec 

(a) Air temperature, 835° Rj air velocity, 82 inch es per second; initial 

drop diameter, 550 microns. 


Figure 3*- Radius, temperature, and mass histories of various fuels. 
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(t>) Air temperature, 1360° R; air velocity, 102 inches per second; 
initial drop diameter, 500 microns. 


Figure 3*- Concluded. 









(a) n-hexane. 


Air temperature,. °R 


Drop diameter. 1.600 microns: air velocity. QO inches per second 


Figure 4.- Experimental vaporization time as a function of air temperature. 








(t) n-decane. Drop diameter, 1,600 microns; air velocity, 90 inches per second. 

Figure 4, - Continued. 
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Normal boiling tamper B.ture, °R 


(a) Air temperature, 900° R; drop diameter, 1,600 to 1,700 microns; velocity 

80 inches per second. 

Figure Experimental vaporization time as a function of normal ho ning 

temperature for various fuels. 
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(b) Air tenperature, 1,450° R; drop diameter, 1,50C 

80 Inches per second. 


Figure 5-- Concluded. 
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(c) n-hexadecane. Air temperature, 1,350° E; air velocity, 

80 inches per second. 


Figure 6.- Concluded. 




Steady- atat© temperature. 


n-Hexane 616 

n-Keptane 069 

Iso-octane 670 

lpoo -n-Octane 7 ia 

n-Nonane 763 

n-Decane 90 s 

n-Dodecane sal 

n-Ietradeeane 950 

n-Hexadecane lpo 7 

<y 900 - n-Octadecane ip 65 


Air ten 5 >«rature - 


Boiling 

tras>erature, °R 
616 

669 

670 

718 

763 

805 

881 


Ingobo (ref. 12) 


n-Octade cane 


n-Hexaaecane 


m3aij 




n-Ietradecane 








n-Dodecane 


n-Decane 




n-Oetane 




Iso-octane 


KTifwawwl 


— — ■ ■ I I P— I— — I 


TmmmMMmirjrMfMMMMmm »»w< 


Air temperature, °R 




MSOT ° 7 ” ta^peratura aa a function of air 

.temperature for various fuels. 
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Batlo of experimental to calculated time for 
20-percent mass VEiporliation 



(a) Ratio for 20-percent mass vaporization. 

Figure 11.- Ratio of experimental to calculated time to reach a percentage mass vaporization 

as a function of drop size for various fuels. 
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(a) Calculations with and without unsteady state considered. 


Figure 12.- Comparison of ratio of experimental to calculated vaporization times with "boiling 

-t-orrm»r>n+.nr>«i -Pn-p* imrinno -PhrIo. Th’nn a -I rp _ SOH mi ^Tfinfl ; R i T ■f'.RTrmPTRtilT'fi . 1 . ^00 ^ . 
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Horaial ‘boiling temperature, °R 

(b) Calculation performed for "zero" thermal conductivity or "peeling 

Figure 12.- Concluded. 
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700 800 900 1POO 

Horraal boiling point, °R 


Figure 15.- Comparison of the physical ignition delay of normal paraffin 
fuels. Air temperature, 1200° R; air velocity, 100 inches per second 
initial diameter, 500 microns; initial fuel temperature, 500° R. 
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Fuel-all ratio, by iwigrt 

Fiprure l6. - Overall fuel-air ratio versus fi n al, adiabatic saturation temperature for^n-heptane 
and n-bexadecane . Air temperature, 1,200' J R; initial fuel temperature } pOu w it. 
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(to) n -hexaclecane . 
Figure 17-- Concluded. 
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(a) Vapor pressure against distance 
from, drop surface at various 
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(b) Temperature distribution in film at 
various times in history. 


(d) Liquid temperature 
against time. 9^ 



(e) Temperature of stoichiometric 

mixture against distance in film. 


particular value of 
phys ical-ignit ion- 
delay time. 


Figure 18.- Development of self -ignition mixture. 
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Figure 19 .- Stoichometric-mixture temperature against time for n-decane; 
self -ignition temperature, 909° R. Air temperature, 1,200°; air 
velocity, 100 inches per second; initial drop radius, 500 microns. 
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